Electron Transport through a Molecular Conductor with Center-of-Mass Motion 
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The linear conductance of a molecular conductor oscillating between two metallic leads is investigated nu- 
merically both for Hubbard interacting and noninteracting electrons. The molecule-leads tunneling barriers 
depend on the molecule displacement from its equilibrium position. The results present an interesting inter- 
ference which leads to a conductance dip at the electron-hole symmetry point, that could be experimentally 
observable. It is shown that this dip is caused by the destructive interference between the purely electronic and 
phonon-assisted tunneling channels, which are found to carry opposite phases. When an internal vibrational 
mode is also active, the electron-hole symmetry is broken but a Fano-like interference is still observed. 
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Molecular electronics has received much attention in the 
past decade, particularly since it became possible to fabricate 
devices in which the active element is a single molecule. 1 - 2 
A fundamental property of molecular conductors is their dis- 
crete electronic spectrum. Although the weak coupling of the 
molecule to the two metallic electrodes leads to the broad- 
ening of the molecular energy levels, their discrete nature is 
maintained. Due to the small size of these molecules, elec- 
tronic correlations are dominant and they lead to interesting 
many-body effects, such as the Coulomb blockade and Kondo 
resonance 3 . These effects have been observed experimentally 
in molecular conductors 4 and other nanostructures. 5 Another 
interesting property of molecules is their flexible nature. They 
have an intrinsic spectrum of internal vibrational modes and, 
when coupled to the electrodes, some molecules acquire ex- 
ternal vibrational modes as well. The excitation of one or 
more of these modes leads to the modulation of the electronic 
energy levels and tunneling barriers between the molecule 
and the electrodes or between different parts of the molecule, 
thus changing the molecular transport properties. These vi- 
brational effects have been observed in a number of recent 
experiments, 6 7 and have been the subject of considerable the- 
oretical investigation 8 . 

In this work, we study the zero bias conductance of a 
molecular conductor model with one relevant electronic en- 
ergy level, both with interacting and noninteracting electrons. 
The molecule is allowed to oscillate between the two elec- 
trodes. This center-of-mass (CM) vibrational mode is treated 
quantum mechanically and leads to an asymmetric modula- 
tion of the tunneling barriers molecule-electrodes. The vibra- 
tional excitation is also coupled to the charge on the molecule. 
This is due to the fact that the chemical bonds inside the 
molecule and the molecule-electrode bonds depend in gen- 
eral on the molecule's charge. The results show an interest- 
ing and unexpected conductance cancellation when an odd 
number of electrons occupy the molecule. It is discussed be- 
low that this cancellation is due to the destructive interference 
between the purely electronic and phonon-assisted tunneling 
channels, which are found to carry opposite phases. 910 In this 
case both channels are elastic. The phonons are virtual, not 
thermal. 
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FIG. 1: A schematic of the system studied in this paper. The 
molecule can oscillate between the two leads around the equilibrium 
position lo, thus modulating the molecule-leads tunneling barriers. 



Figure ^ schematically depicts the system analyzed in this 
work. The molecule can oscillate between the source and 
drain electrodes, thus modulating the tunneling barriers. In 
our calculations, this modulation and the electron-vibration 
coupling are expanded up to the linear term. 1112 The elec- 
tronic part of the system is modelled using the Anderson im- 
purity Hamiltonian. The total Hamiltonian can be written as 
H = Hm + ^lcads + -f*M-ieads, where Hm is the Hamiltonian 
of the molecule, 
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The first term represents the energy of the relevant molecular 
orbital controlled by the gate voltage, the second term repre- 
sents the Coulomb repulsion between the electrons occupying 
the molecular orbital, the third term couples the vibrational 
excitation to the net charge on the molecule (aJ and a are the 
phonon creation and annihilation operators), and the fourth 
term represents the vibrational energy. Hi ca ds describes the 
two leads modeled here as semi-infinite ideal chains, 
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where c\ i(y {c\, i(T ) creates an electron with spin a at site i in the 
left (right) lead, t is the hopping amplitude in the leads and 
the energy scale (t = 1). flM-leads connects the molecule to 
the leads, 
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FIG. 2: (a) NRG and (b) ED+DE results for G as a function of V g 
in the Kondo Regime for a — and increasing the electron-phonon 
coupling strength (X/u>o = 0.0, 0.4, 0.8, and 1.2). For small X/ujo, 
the standard Kondo resonance is obtained with renormalized param- 
eters. In the limit of large X/ujq, the charge Kondo effect is obtained. 
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where d\ creates an electron with spin a in the molecule, t' 
is the hopping parameter between the molecule and the first 
site of each lead, and a is a parameter that carries the depen- 
dence of t' on the molecule displacement from its equilibrium 
position x (note the opposite signs in this dependence for the 
two leads). This displacement can be written in terms of the 
phonon operators as x = (a+a^). The total Hamiltonian is in- 
variant under the particle-hole and (a — > —a) transformation. 
In the results shown, unless otherwise stated the following set 
of parameters was used (U = 1.0, t' — 0.2, luq — 0.2) while 
A and a were variedtAiThe value of u>q was fixed since its in- 
crease or decrease would simply produce the opposite effect 
of increasing or decreasing A and/or a. 

Using the Keldysh formalism, 14 the zero bias and 
zero temperature conductance can be written as G = 
^-|i 2 Gz r (ej?)| 2 [p(ei?)] 2 , where G; r is the Green's function 
that propagates an electron from the left to the right lead and 
p(cf) is the density of states in the leads at the Fermi level. 
Note that at zero bias only elastic processes can be observed. 
The Green's functions are calculated using exact diagonaliza- 
tion supplemented by a Dyson equation embedding procedure 
(ED+DE). (See Ref.| 15] for a full description of the method). 

It is useful to start by studying briefly the model when the 
tunnel barriers do not depend on the vibrational excitation 
{a = 0). This case was previously studied using numeri- 
cal renormalization group (NRG) techniques. 16 In this case, 
the model can be mapped into an effective electronic Hamil- 
tonian with renormalized parameters (U = U — 2A 2 /wo, 
V g = V g + X 2 /u Q and ? oc t' exp(^)) 16 - 11 . FiguresEJi and 
b show the results in the interacting electrons case obtained 
using ED+DE and NRG, respectively. For weak electron- 
phonon coupling, the standard Kondo resonance is obtained 
with renormalized parameters. In the strong electron-phonon 
coupling limit, the 'charge Kondo effect' is obtained. Note 
that the results obtained using ED+DE are very similar to the 
NRG results, clearly capturing the essence of the problem. 
Since NRG methods cannot be applied to the CM oscillations 
studied here, and ED+DE appears equally accurate, the results 
presented below were obtained using this last technique. 




FIG. 3: (a) G as a function of V g in the Kondo regime for a = 
(solid line) and a = 0.4 (dashed line), A = 0.2 in both cases. In the 
first case, the usual Kondo peak with reduced width is obtained. In 
the second case, a conductance dip is obtained, (b) Average occupa- 
tion of the molecular level for the same set of parameters. Note that 
the charging behavior is almost the same in both cases, (c) G vs. V g 
in the absence of Coulomb repulsion (U = 0), for a = (solid line) 
and a = 0.4 (dashed line), A = 0.1 in both cases. The conductance 
dip is also obtained in this case, thus the physical mechanism behind 
this effect does not depend on the electron-electron interaction. 



Figure |3 contains the main results of this work. Figure [3^ 
shows the conductance in the interacting electrons case. For 
a = 0, the conductance simply shows a Kondo resonance 
peak with reduced width. However, when a ^ 0, a conduc- 
tance dip is obtained when an odd number of electrons occupy 
the molecule. Figure[3p provides the average occupation (rid) 
of the molecular orbital where it can be clearly seen that the 
charging behavior is almost the same in both cases. Note that 
for the usual Friedel sum rule 18 G = sin (^(n^)) 

is not satisfied and this can be an indication of a non-Fermi liq- 
uid behavior. Figure^ shows the conductance in the absence 
of Coulomb repulsion (U — 0): the same effect is obtained 
as for the case a ^ 0. The conductance cancellation does 
not depend on the electron-electron interaction, which agrees 
with the qualitative explanation presented below. 

Figure0]shows the conductance as a function of V g for dif- 
ferent values of A. The dip becoming more pronounced as A 
increases i.e. as the average number of phonons in the ground 
state increases. The inset shows the results obtained by trun- 
cating the phonon Hilbert space at different maximum number 
of phonons (-/V p h). In all the calculations, iV p h = 7 was used 
unless stated otherwise. The qualitative effect of conductance 
cancellation is obtained all the way down to iV p h = 1, al- 
lowing us to study larger clusters and reduce size effects to 
intuitively understand the origin of the dip. 

In figure |5] an explanation of the conductance dip is pre- 
sented. The reasoning starts by noting that i?M-icads in 
Eq|3] can be rewritten as a sum of two channels contribut- 
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FIG. 4: G as a function of V g in the Kondo regime for a — 0.4 and 
A =0.15, 0.20, and 0.25 (dashed, dotted and solid lines respectively). 
The dip becomes more pronounced as A increases. The inset shows 
the convergence of the results with the maximum number of phonons 
(JVph). Note that the qualitative effect of conductance cancellation is 
preserved all the way down to N p h = 1. 



ing to the overall molecule-leads connection. The first term, 
^ Tier {dt c io<j + dJ-CrOer + h.c), represents the purely elec- 
tronic tunneling between the molecule and the two electrodes. 
The second term, t'a{a + a { )^ rJ {d\c rQa — d\cio a + h.c.) , 
represents a phonon assisted tunneling channel, i.e. the elec- 
tron absorbs (emits) a phonon upon entering the molecule and, 
then, emits (absorbs) a phonon upon leaving. Note that both 
channels are coherent and elastic. The number of phonons in 
the system does not change. Fig. |5t shows a schematic of 
the two channels. These channels were studied separately by 
keeping only the relevant term in Hm- leads- The conductance 
and the phase carried by each channel were calculated. Fig. 
13) shows the conductance of the separate channels. Fig. |5J; 
contains the conductance when both channels are active i.e. 
when both terms are included in IfM-icads- Fig- 01 shows 
the phase difference between the two channels. Note that for 
V g = —U/2, the conductance of each of the channels is 2e 2 /h 
and the phase difference is 7r, leading to a perfect cancellation 
in the overall conductance. 10 This interference effect is inde- 
pendent of the electron-electron interaction and, thus, the can- 
cellation should still be present for U = as already shown. 
The dip becomes more pronounced as A increases, increasing 
the average number of phonons in the ground state. 

The stability of the dip is tested by adding an internal vibra- 
tional mode which leads to the symmetric modulation of the 
tunnel barriers to the leads (breathing mode). To account for 
this mode, the following terms were added to the Hamiltonian 

H' = A'(l -n d )(b + tf) +uj' tfb + 

t'a'(b + 6 f ) {dUwa + dic r0a + h.c), (4) 

(7 

where the first term represents the electron-phonon coupling, 
the second term represents the breathing vibrational energy, 
and the third term represents the subsequent modulation of the 
tunnel barriers. The results are shown in Fig[6] When only the 
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FIG. 5: (a) A schematic representation of the two conductance chan- 
nels, the purely electronic tunneling represented by the two upper 
arrows and the "phonon-assisted tunneling" channel represented by 
the two lower arrows, (b) Partial conductance when only one of 
the channels is active. The dotted line shows the conductance of 
the purely electronic tunneling channel Ge, while the dashed line 
shows the conductance of the "phonon-assisted tunneling" channel 
G\. (c) Conductance when both channels are active, (d) Phase dif- 
ference A$ of the two channels. Note that A$ for all values of V g is 
close to 7r thus leading to destructive interference. In particular, for 
V g = -U/2, A$ = 7T and G E = Gi = 2e 2 /h, thus leading to a 
perfect cancellation in the overall conductance. (A = 0.2 , a — 0.4). 
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FIG. 6: The solid line shows the conductance of the molecule when 
a breathing vibrational mode is active (no CM motion). The particle- 
hole symmetry is broken as expected and no conductance dip is ob- 
tained. The dashed line shows the conductance when both breathing 
and CM vibrational modes are active. The combined effects of the 
two modes lead to a Fano-like interference. The breathing mode pa- 
rameters used are A' = 0.2, u>' = 0.3 and a' — 0.3. 



internal mode is active (solid line), the electron-hole symme- 
try is broken but no dip is observed. This agrees with previous 
results 19 obtained using NRG calculations. In the case where 
both vibrational modes are active (dashed line), the dip ap- 
pears. The combined effect of conductance cancellation and 
electron-hole asymmetry leads to a Fano-like interference. 

The finite-size effects on the results are shown in Fig. 01 
where the convergence of the conductance with the size of the 
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FIG. 7: (a) Convergence of the conductance with the size of the 
exactly-solved cluster L. The solid line, dashed line, and the circles 
show the results obtained using L = 3, 7, and 11 respectively. In the 
three cases, a maximum number of three phonons was used, (b) The 
real and (c) imaginary parts of the isolated cluster Green's function 
giz that propagates an electron from the left to the right ends of the 
cluster (L — 3 and V g = —{7/2). Note that both parts are equal 
to zero at the Fermi level (located at uj — 0.0). Thus, the origin of 
the conductance dip can be traced back to the exactly-solved cluster. 
A = 0.2 and a = 0.4 in the three figures. 



exactly-solved cluster is presented. Note that increasing the 
cluster does not change the qualitative effect of the conduc- 
tance dip. Moreover, the origin of the conductance dip can be 



traced back to the exactly-solved cluster by studying Green's 
functions before the embedding process. Figures [7p and c 
show the real and imaginary parts of the Green's function g\ T 
that propagates an electron from the left to the right end of the 
cluster for a ^ and V g = —{7/2. Both parts are zero at the 
Fermi level (u> = 0). For a — (not shown here), g\ r has a 
pole at the Fermi level and the system is perfectly conducting. 
When a is turned on, the pole splits into two, one below and 
one above the Fermi level thus causing the zero conductance. 

In conclusion, the zero temperature electron transport 
through a molecular conductor with center-of-mass motion 
was studied numerically for interacting and noninteracting 
electrons. The results present an interesting conductance dip 
when an odd number of electrons occupy the molecule. It is 
argued that this dip is caused by the destructive interference 
between the purely electronic and phonon-assisted tunneling 
channels, which were found to carry opposite phases. When 
an internal vibrational mode is also active, the particle-hole 
symmetry is broken but a Fano-like interference is still ob- 
tained. The conductance cancellation would best be observed 
on a broad conductance peak such as the Kondo peak which 
is broader than the resonant tunneling peak. 
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